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The effect of some possible non standard WIMP velocity distributions, like the Debris Flows 
recently proposed, on the direct dark matter detection rates is investigated. We find that such 
distributions may be deciphered from the data, especially if the time variation of the event rates 
due to the annual motion of the Earth is observed. 
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I. INTRODUCTION 



n b = 0.0456 ± 0.0015, fi CDM = 0.228 ± 0.013, fl A = 0.726 ± 0.015. 
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■ The combined MAXIMA-1 Q, BOOMERANG 0, DASI and COBE/DMR Cosmic Microwave Background 
| (CMB) observations [H imply that the Universe is flat [|| and that most of the matter in the Universe is Dark [1, @], 
i.e. exotic. Combining the data of these quite precise experiments one finds: 

O [ Since any " invisible" non exotic component cannot possibly exceed 40% of the above r2 CDM Q , exotic (non baryonic) 
matter is required and there is room for cold dark matter candidates or WIMPs (Weakly Interacting Massive Particles) . 

Even though there exists firm indirect evidence for a halo of dark matter in galaxies from the observed rotational 
curves, see e.g the review 0, it is essential to directly detect such matter. The possibility of such detection, however, 
t-H ■ depends on the nature of the dark matter constituents and their interactions. 

| Since the WIMP's are expected to be extremely non relativistic, with average kinetic energy (T) sa 
50 keV(m wimp/100 GeV), they are not likely to excite the nucleus. So they can be directly detected mainly via 
the recoiling of a nucleus (A,Z) in elastic scattering. The event rate for such a process can be computed from the 
following ingredients [lOj : i) The elementary nucleon cross section. This most important parameter will not, however, 
be the subject of the present work. We will adopt the view that it can be extracted from the data of event rates, if 
(ST and when such data become available. From limits on the event rates, one can obtain exclusion plots on the nucleon 
cross sections as functions of the WIMP mass, ii) knowledge of the relevant nuclear matrix elements obtained with 
as reliable as possible many body nuclear wave functions, iii) knowledge of the WIMP density in our vicinity and its 
velocity distribution. 

In the standard nuclear recoil experiments, first proposed more than 30 years ago (Tlj . one has to face the problem 
that the reaction of interest does not have a characteristic feature to distinguish it from the background. So for the 
expected low counting rates the background is a formidable problem. Some special features of the WIMP-nuclear 
interaction can be exploited to reduce the background problems. Such are: 

i) the modulation effect: this yields a periodic signal due to the motion of the earth around the sun. Unfortunately 
this effect, also proposed a long time ago [12] and subsequently studied by many authors [ToL [l3Tj2l| . is small and 
becomes even smaller than 2% due to cancelations arising from nuclear physics effects, 

ii) backward-forward asymmetry expected in directional experiments, i.e. experiments in which the direction of the 
recoiling nucleus is also observed. Such an asymmetry has also been predicted a long time ago [22j |. but it has not 
been exploited, since such experiments have been considered very difficult to perform. Some progress has, however, 
has recently been made in this direction and they now appear feasible [2214341 ]. In such experiments the event rate 
depends on the direction of observation. In the most favorable direction, opposite to the sun's direction of motion, 
is comparable to the standard event rate. The sensitivity of these experiments for various halo models has also been 
discussed [27],[28j]. Furthermore we should mention that in such experiments [23], [34| all events are counted. If some 
interesting events can be found, they can be established by further analyzing them by the direction of the observed 
recoils. 

An essential ingredient in direct WIMP detection is the WIMP density in our vicinity and, especially, the WIMP 
velocity distribution. The dark matter in the solar neighborhood is commonly assumed to be smoothly distributed 
in space and to have a Maxwellian velocity distribution. Some of the calculations have considered various forms of 
phenomenological non symmetric velocity distributions f35T - l38j [l9l . [26l . [27j and some of them even more exotic dark 
matter flows like the late infall of dark matter into the galaxy, i.e caustic rings [39M43l | and Sagittarius dark matter 
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In addition to the above models very recently it was found that the velocity distributions measured in high res- 
olution numerical simulations exhibit deviations from the standard Maxwell-Boltzmann assumption, especially at 
large velocities [45l 14(1] . Furthermore a distinction was between a velocity structure that is spatially localized, such as 
streams (47i |48|. and that which is spatially homogenized, which was designated as "debris flow" [49j | . Both streams 
and debris flows arise from the disruption of satellites that fall into the Milky Way, but differ in the relative amount 
of phase-mixing that they have undergone. Implications of streams (50j and the debris flows in direct dark matter 
searches have also been considered (51). 

In the present paper we will discuss in some detail the effect of these debris flows on the event rates of direct dark 
matter experiments for a variety of targets such as those employed in XENON10 (52|,CoGENT[53|, DAMA(54l. |55| 
and PICASSO (5(| We will also study the effect of these flows on the time variation of these rates due to the 
annual motion of the Earth 58] (modulation effect) as a function of the energy transfer and the WIMP mass and 
compare them with the standard M-B distribution. The effects of debris flows in directional experiments will be 
studied elsewhere. 



II. THE FORMALISM FOR THE WIMP-NUCLEUS DIFFERENTIAL EVENT RATE 

The most interesting quantity which depends on the velocity distribution is the quantity g(v m in)- For the M-B 
distribution in the local frame it is defined as follows: 

g(v min , v E {a)) = 1 3 f™ e-^+^^+^y^vdvdn (1) 

For isotropic debris flows [5l[ it is given by: 

( ( \\ f ft \ r 2v fl 7^7aj > v fiow-v E {a) <v <v flow +v E {a) 

9(v min ,v E (a)) = jf^ —dv, f(v) = { otherw . se (2) 

These functions are shown in Fig. [TJ 

In what follows we will find it useful to useful to expand g(v m i n ,VE(a:)) in in powers of 5, the ratio of the Earth's 
velocity around the sun divided by the velocity vq of the sun around the galaxy (220km/s). Keeping terms up to 
linear in 6 ss 0.135 and expressing everything in dimensionless variables we obtain: 

v g(v min ,v E {a)) =ty (x) + ty 1 (x)cosa, x =^^- (3) 

vo 
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where "Jo (2; represents the quantity relevant for the average rate , ^i{x, which is proportional to 5, represents the 
modulation and a is the phase of the Earth (a = around June 3nd) . In the case of the flows they were derived from 
the semi-analytic approximations of simulations as discussed by Spergel and co-workers [5lj . 
In the case of a M-B distribution take the following form: 

*o(a?) = ^(erf(l - x) + erf(x + 1) + crfc(l - y yeec ) + erfc(t/ ye8C + 1) - 2) (4) 

1 /— erffl — x) — evilx + 1) — erfcfl — y OBC ) — erfc(?/ csc + 1) 
= 2% 2 

e -(^-i) 2 e -(^+i) 2 e -feo S c-i) 2 e -(yo S o+i) 2 



1 (5) 
V" V" V" V" / 

where erf(x) and erfc(x) are the error function and its complement respectively y esc w 2.8uo is the escape velocity. 
For isotropic debris flows one finds: 



*"«•'■' = ^ Vf-Kx<l + Vf , Vf = VJ ^ (6) 

x > 1 + yf 



< x < y f - 1 




< x < y/ - 1 

; — ?y/ -l<x<l- ?y/ , 1// = ^ (7) 
a; > 1 + yf ' " 

We note that the variable a; depends on the nuclear recoil energy Er as well as the WIMP-nucleus reduced mass. As 
we shall see below there is an additional dependence of the rates on Er coming from the nuclear form factor. 

At Earth-frame velocities greater than 450 km/s, debris flow comprises more than half of the dark matter at the 
Sun's location, and up to 80% at even higher velocities [5l[. In the VL2 simulation, the combination of debris flows 
and standard M-B is very well fit by the function 

£(I ) = 0.22 + 0.34 ( etf (4g-i§) +1 ) (8) 

This function is exhibited in Fig. [2] In this case we find: 

*i{x) -> (l-e(x))V? IB (x)+e{x)y{(x), i = 0, 1 (9) 

The functions ^o(x) and 1 $i(x) are exhibited in Fig. [3J As expected in the case of the flows ^q(x) falls off linearly 
for large values of x. Note that in all cases ^i(x) takes both positive and negative values, which affects the location 
of the maximum of the modulated rate as a function of a, depending on the target and the WIMP mass. We will 
explore this effect of the different distributions in direct experiments searching any time dependence of the rates. 

Once these functions are known the formalism to obtain the direct detection rates is fairly well known (see e.g. the 
recent reviews [H, So we will briefly discuss its essential elements here. The differential event rate can be cast 

in the form: 

dR dRo dH 

^ U = ^ U + ^ UC0Sa (10) 

where the first term represents the time averaged (non modulated) differential event rate, while the second gives the 
time dependent (modulated) one due to the motion of the Earth (see below) . Furthermore 

1 dt 



dRo , 


Px 


m t . 


( ^ 


dE R lA ~ 


m x 


Am p n 


KPp 


dH 


Px 


m t 


( El 


dE R U " 


m x 


Am p " 


Kp p 



V< v 2 >A 2 



Qo(A)du' 



2 VZv^A 2 -^-^ (11) 
Qo{A) du 



4 



0.8 
0.7 
0.6 
0.5 
0.4 




FIG. 3: The functions ^o(x) and vE'i(x) as a function of x — v m i n /vo. Note that these functions have been 
computed at a — 0, i.e. the average local velocity. Note also that the variable x depends on the nuclear recoil energy 
En as well as the WIMP-nucleus reduced mass. Otherwise the labeling of the curves is the same as that of Fig. Q] 



with with \i r (fi p ) the WIMP-nucleus (nucleon) reduced mass, A is the nuclear mass number and a n is the elementary 
WIMP-nucleon cross section. m x is the WIMP mass and mt the mass of the target. Furthermore one can show that 



(12) 



with a = (^/2fj, r bvo) 1 , vq the velocity of the sun around the center of the galaxy and b the nuclear harmonic oscillator 
size parameter characterizing the nuclear wave function, u is the energy transfer Q in dimensionless units given by 

E " , Q (A) = [mpAb 2 ]- 1 = 40^- 4 / 3 MeV (13) 



Qo(A) 

and F(u) is the nuclear form factor. In the present calculation they were obtained in context of the nuclear shell 
model in the spirit of[6l| (for the spin induced process see, e.g. [HJ, The form factor is important in the case of 
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a heavy target and large WIMP mass, i.e. for large recoil energies (see Fig. 0]). 

Note that the parameter a depends both on the WIMP , the target and the velocity distribution. Note also that 
for a given energy transfer Er the quantity u depends on A. 
Sometimes one writes the differential rate as: 

dR p x m t ( fi r \ 2 r- 2 - 2 1 / dt 



dEft m x Am 



-a n (g) VZ^A*^ (£(l + i7 (a ^7Q^4))cos«) (14) 



In this formulation H(ay / En/Qo(A)), the ratio of the modulated to the non modulated differential rate, gives the 
relative differential modulation amplitude. It coincides with the ratio 4 , i(ay / E^/Qo(A)) /^o(a^/ Er/Qo(A)), i.e. it is 
independent of the nuclear form factor and depends only on the reduced mass and the velocity distribution. It is thus 
the same for both the coherent and the spin mode. Note that it can take both positive and negative values, which 
affects the location of the maximum of the modulated rate as a function of a. For the convenience of the analysis of 
experiments, however, we will present our results in the form of Eq. 1111 



III. SOME RESULTS ON DIFFERENTIAL RATES 



We will apply the above formalism in the case of I and Na, which are components of the target NAI used in the 
DAMA experiment [HI, HI] and Ge employed, e. g, b y the CoGeNT experiment The results for the Xe target [52[ 
are similar to those for I and for the 19 F target [56l are similar to those for Na . The differential rates %q\a an d 

f° r each component (A = 127 and A = 23) and for A = 73 are exhibited in Fig. ISlfTUl The nuclear form factor 
has been included (for a heavy target, like 127 I or 131 Xe, its effect is sizable even for an energy transfer[58| of 10 keV, 
see Fig. 

The introduction of debris flows makes a small contribution at low energy transfers. As expected [5l| it tends to 
increase the differential rate at high energy transfers. This is particularly true for light small WIMP-nucleus reduced 
mass (see Figs and [9]). One, however, does not see any particular signature in the shape of the resulting curve. 
Furthermore the event rate in this region is about five times smaller than the maximum. One, however, observes an 
interesting pattern concerning the time varying (modulated) part of the rate (see Figs lr?l 151 and fTT])) . For a heavy target, 
like 127 I or 131 Xe, it is not surprising that, for WIMPs with relatively large mass, the modulation becomes negative, 
i.e. the rate becomes minimum in June 3nd, for all models considered here. For low WIMP masses, however, the sign 
of the modulation due to the flows is opposite to that of the M-B distribution. Thus the use of the light target nucleus 
19 F, combined with the low detection threshold of 1.7 keV for recoil nuclei, makes PICASSO particularly sensitive to 
low mass dark matter particles and give s it also some leverage in the low mass region of the spin independent sector. 
The present stage of the experiment [631] is a ppr oaching the sensitivity to challenge or confirm the claims of seasonal 
modulations by the DAMA[55| and CoGeNT[53j experiments. A similar situation arises in the case of an intermediate 
target, like 73 Ge. Here the M-B distribution yields a negative value only for very low energy transfers. The situation 
becomes most interesting in the case of a light target, see Fig. [5J Here, with the possible exception of quite low 
energy transfers, which perhaps are below or very near threshold, the M-B distribution yields a positive modulation 
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amplitude, i.e. a maximum on June 3nd, while the result of debris flows is to cause a change in sign as one moves to 
high energy transfers. Also in this case the modulation amplitude tends to increase as the energy transfer increase, 
while the corresponding contribution due to the M-B distribution tends to decrease. We should remark though that 
the total rate (average+modulated) tends to decrease at high momentum transfers. We should also stress that we 
have presented here the absolute modulate rate (events per kg target per year). The relative modulated amplitude 
(the ratio of the time varying rate divided by the time averaged) maybe larger. 

The above results, as we will see in the next section, have important implications in the total event rates. 

Sometimes, as is the case for the DAMA experiment, the target has many components. In such cases the above 
formalism can be applied as follows: 

——1^4 — > -^37) \ A n u Ui > = ^ ne f rac tion of the component A\ in the target (15) 
We will not, however, pursue such an analysis. 



IV. SOME RESULTS ON TOTAL RATES 



For completeness and comparison we will briefly present our results on the total rates. Integrating the differential 
rates discussed in the previous section we obtain the total rate R, adding the corresponding time averaged rate Rq 
and the total modulated rate H, given by: 

R = R + H = (BL \ V< v 2 >A 2 a n t (1 + h cos a), (16) 

m x Am p \fi p J 

with 

t = / — du, h = - I —du. (17) 

jQ th /Qa(A) du t J Q th /Q Q (A) du 
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(c) (d) 
Q -^keV 

FIG. 6: The differential rate asa function of the recoil energy for a heavy target, e.g. 127 I assuming a nucleon 
cross section of 10~ 7 pb. Panels (a) (b), (c) and (d) correspond to to 5, 20, 50 and 100 GeV WIMP masses. 

Otherwise the notation is the same as that of Fig. [1] 



y max is the maximum velocity allowed by the distribution and Qth{A) is the energy cut off imposed by the detector. 

The obtained results for quantities Ro and h are exhibited in Figs [TT1fT51 assuming a nucleon cross section of 10 _7 pb. 
In the case of a heavy target the average event rate attains the maximum value of 30 events per kg of target per 
year at a WIMP mass of 25 GeV, while for heavy WIMPS it eventually falls to about 5 kg/y at 500 GeV(to a good 
approximation it falls inversely proportional to the WIMP mass above the 200 GeV). For an intermediate target we 
get 15 kg/y at 25 GeV, with an asymptotic value of 4 kg/y. For a light target the maximum becomes 2.5 kg/y at 20 
GeV. Again the asymptotic value at 500 GeV is about 1/5 of the maximum. 

It is clear that, as far as the time average rates Rq are concerned ,the debris flows do not exhibit any characteristic 
signature to differentiate them from the standard M-B distribution. The relative modulation amplitude h, however, 
exhibits a very interesting feature, namely, if caused by the flows, it is negative for all targets, even for the light ones, 
and in the entire WIMP mass range (minimum in June). On the other hand if it is caused by the M-B distribution it 
is positive in the case of light targets regardless of the WIMP mass. It is also positive for intermediate/heavy targets, 
if the WIMPs are relatively light. Then the maximum occurs on June 3nd as expected. It becomes negative only for 
relatively heavy WIMPs. Thus it is an experimental challenge to measure the small time dependence of the event 
rate with a relative difference between the maximum and the minimum of 2h ss 4%. From such data on both light 
and heavy targets, if and when they become available, one may may be able: i)to get a hint about the size of the 
WIMP mass and ii) infer the existence of flows. 

V. DISCUSSION 

In the present paper we first obtained results on the differential event rates, both modulated and time averaged, 
focusing our attention on the effects of debris flows. We found that: 

• The flows indeed enhance the time averaged rates at relatively high energy transfers compared to the M-B 
distribution. All rates, however, fall as the energy transfer increases. This fall is only partially due to the 
velocity distribution. It is also caused by the nuclear form factor, in particular in the case of heavy targets. 
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(c) (d) 
Q -^keV 

FIG. 7: The differential rate as a function of the recoil energy for a light target, e.g. 23 Na assuming a nucleon 
cross section of 10 _7 pb. Panels (a) (b), (c) and (d) correspond to to 5, 20, 50 and 100 GeV WIMP masses. 

Otherwise the notation is the same as that of Fig. [T] 



• In view of the dependence of above rates on the unknown WIMP mass, one does not have a clear signature to 
differentiate the debris flows from the standard distribution. 

• The differential modulated rates provide such a signature, the sign of the modulation amplitude, which deter- 
mines the position of the maximum. The debris flows tend to favor a negative sign (minimum on June 3nd), 
while the standard WIMPs favor a positive sign when the target is light or even when the target is heavy but 
the WIMP is light (maximum on June 3nd). 

We then proceeded and calculated the total event rates as functions of the WIMP mass. We presented results obtained 
with a zero threshold ener gy. For higher threshold energies we expect the debris flows to be suppressed a little less 
than the standard WIMPs 58], since the differential event rates associated with the former attain smaller values at 
low energy transfers. The time averaged rates are affected by the debris flows, but one does not find a characteristic 
feature to differentiate the debris flows from the standard WIMPs. Such a feature is provided by the sign of the 
modulation amplitude. The debris flows tend to favor minimum in June 3nd (opposite to that expected), while for 
the standard WIMPs it is positive (maximum on June 3nd) in the case of light targets, regardless of the WIMP mass, 
or for heavy targets for relatively light WIMPs. 

In conclusion a measurement of the modulation amplitude expected at the scale of 2h « 4% may shed light on the 
scale of the WIMP mass and determine the presence of any non standard velocity distributions. These issues maybe 
settled better, if data on directional experiments become available. Such theoretical explorations are currently under 
stydy. 

Acknowledgments 

The author is indebted to the PICASSO collaboration for partial support of this work and to Viktor Zacek and Ubi 
Wichoski for their kind hospitality in Montreal and SNOWLAB. This work was also partially supported by UNILHC 
PITN-GA-2009-237920. 



9 




(c) (d) 
Q ->keV 

FIG. 8: The differential rate ^, asa function of the recoil energy for a light target, e.g. 23 Na assuming a nucleon 
cross section of 10~ 7 pb. Panels (a) (b), (c) and (d) correspond to to 5, 20, 50 and 100 GeV WIMP masses. 

Otherwise the notation is the same as that of Fig. Q] 
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FIG. 11: The total rate Rq (top panels) and the relative modulation amplitude h (bottom panels) as a function of 
the WIMP mass in GeV in the case of a heavy target 127 I at zero threshold. Otherwise the notation is the same as 

that of Fig. CD 





FIG. 13: The same as in Fig. [TT]for an intermediate target, e,g. 73 Ge 



